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Abstract 

The influence of the dissipative terms on the conditions of forma- 
tion and the characteristic parameters of shock waves in relativistic 
nuclear collisions is investigated for three types of equation of state 
(non linear QHD-1, resonance gas and lattice QCD). Energy and ve- 
locity proflles are obtained in a one- dimensional model; the duration 
of the shock phase and width of the shock front are calculated. It is 
shown that the presence of a phase transition results in a strong en- 
hancement of the width of the shock front, which results in an increase 
of transparency. This effect, combined with the fact that the nuclei 
have a flnite size, prevents the energy density to rise to its maximum 
value (full stopping) as would be predicted by a non dissipative shock 
model. 

1 Introduction 

Among the possible approaches for the description of relativistic heavy ion 
collisions, hydrodynamics has obtained some success in reproducing experi- 
mental data [1] - [6]. It is known from nonrelativistic calculations [7] - [10] 
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that dissipative effects can have an important inffuence on the observabies. 
On the other hand, reiativistic effects are thought to be essentiai for a good 
description of nuciear coiiisions above f GeV energy. However, a fuiiy reia- 
tivistic hydrodynamic code inciuding dissipative terms yet practicaiiy does 
not exist (see however [f f]). The probiems hampering this deveiopment are 
of two different natures. First, there is a numericai one: The Lorentz invari- 
ance now requires the time and space coordinates to be treated on an equai 
footing. This introduces time derivatives in the hydrodynamical equations 
which cannot anymore be solved by the usual explicit algorithms. Secondly, 
some doubts have been emitted concerning the validity of the Navier-Stokes 
equations for reiativistic problems, as these may exhibit unstable and acausal 
modes (see e.g. [12] - [16]). 

In this paper we would like to make a first step in the study of reiativis- 
tic dissipative hydrodynamics by addressing the question of the structure of 
reiativistic shock waves in the presence of transport terms. This simplified 
picture of the early stage of the collision allows to find an analytical solu- 
tion and permits to gain insight into the basic physical processes at hand. 
The solution of this problem may then be used as a simple picture for the 
early stage of the collision and as input for the subsequent expansion phase. 
The model yields energy density and velocity profiles, the duration of the 
compression phase and an analytical expression for the width of the shock 
front. If the width of the shock front is larger than the size of the nuclei, the 
notion of shock looses its meaning. It therefore serves as a criteria to what 
extent hydrodynamical models are applicable. Moreover, it provides a basic 
check in the development of a more refined hydrodynamical code and a test 
of some general features of the equation of state. 

2 Basic formalism 

2.1 Reiativistic Navier-Stokes equation 

The reiativistic Navier-Stokes equations [17, 18] are derived from the relation 
of conservation of the baryon number and energy-momentum 
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The baryon current and energy momentum tensor take the general form: 

= pu^ (3) 
T^" = eu^u" - P^^" + 2f]a^'' + CIW^^" + q^u" + q'u^ (4) 

where = (7,7^) is the hydrodynamical 4- velocity, A^'' = g^^ — u^u^ , 
g^^ = diagil, —1, —1, —1). rj and C, are the shear and bulk viscosities; is 
the heat current and is given by 

< = Xl^^'d.iT) -Tu^] = -K^^'d^iulT), (5) 

A is the thermal conductivity and K is related to it through 

K = pT^\l{e + P). (6) 

(T^^ is the shear tensor and 6 is the 4-divergence of the velocity; 

e = dX (8) 

The thermodynamical properties (e, P, p) and transport properties (A, 
?7, of the medium may be given by phenomenological models or derived 
from an underlying microscopical theory (see below). 



2.2 Equation of state 

In the following, four cases were considered for the equation of state (eos): 
In order to have a simple model, we first assumed an equation of the form 
P = (?e with a constant sound velocity, t.g. 1/3. 

We also took the equation from the Walecka model (here with the ex- 
change of uj and a mesons only) [27]. We used the non linear version of the 
QHD-l model, which gives values for the effective mass and incompressibility 
coefficients compatible with the experimental values. The parameters 

m„ = 550MeV; = 783MeV; g^ = 7.947; g^ = 6.706; 

B = -0.17788 10"^ C = 0.39674 10"^ (9) 

yield [28] 

m* = 0.85; K^o = 210MeV, (10) 
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Figure 1: Speed of sound (? for the non-linear QHD-1 model 



in agreement with the current experimental values. 

The use of the non-linear QHD-1 model has the advantage that one can 
derive the relativistic Navier-Stokes equations, the eos and the transport 
coefficients in a consistent manner (see below). 

Next, we used the resonance gas model [29]. In the version calculated 
here, all particles and resonances of the Particle Data Book are taken to 
contribute to the pressure and energy density. 
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where i is running over the particle and resonance states, Ei is the energy 



Figure 2: Speed of sound (? for the bootstrap equation of state 



of particle z, is the degeneracy factor, Bi the baryon number and Si the 
strangeness, and /i^ are the chemical potentials for baryon number and 
strangeness respectively, and p\, and are the corresponding densities. 

The conservation of strangeness is built in by introducing a strange chem- 
ical potential p^ which is adjusted in such a way that the strangeness of the 
system is vanishing. This strange chemical potential grows with baryonic 
density and temperature (c/Fig 3). 

It is interesting to note that, at high baryon densities, p^ exceeds the kaon 
mass. This leads to a condensation of the kaons which we take into account 
by adding a condensate contribution nj^", which is given by 

= -Psil^b, T, fis = mK)0{fis - mx) (14) 

This contribution strongly increases at high temperature and baryonic den- 
sities (cf Fig 4) 

Finally, we considered an eos derived from lattice QCD calculations in- 
cluding a first order phase transition from a quark-gluon plasma towards 
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Figure 3: strange chem- 
ical potential for boot- 
strap eos 



nuclear matter. The lattice data are reproduced by a phenomenological fit- 
ting formula [26]: 

P = cl{e)e 

4(8) = [a + /3tanh (7 In (^) + S)] [l - ~ ] (15) 



Two sets of parameters were used (see Table) 
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/3 
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set (1) 


5/21 


2/21 


0.24 


1.05 


0.3 


0.73 


3.0 


set (2) 


4/15 


1/15 


0.24 


1.05 


0.25 


0.73 


1.5 



In Figs. 1, 2, 5, we show the behavior of = P/e as a function of the 
temperature and of the chemical potential for the various equations of state 
considered. As we will see later, is the shock velocity in the ultrarelativistic 
limit. 
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Figure 4: density of 
kaon condensate for 
bootstrap eos 
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2.3 Transport properties 

The volume viscosity ( is in general negligible, and it plays in the Id case a 
role similar to that of the shear viscosity (c/eqs. (4-7)), so that it will not 
be considered in the following. We will consider three cases for the thermal 
conductivity A and the shear viscosity rj. 

In order to keep things simple, we first chose constant transport coeffi- 
cients. 



Second, in the Walecka model we used a fit of the thermal conductivity 
and shear viscosity as provided by [31]. These transport coefficients can also 
be calculated from the Walecka a - io model in a consistent manner with eos 
[33], which gives values somewhat larger but qualitatively similar to those of 
[31]. In Fig. 6, we show a plot of the viscosity ?y, the thermal conductivity A 
as well as of the parameter K (Fqn. 6 ) entering the equation for the shock 
profile as coefficient of dz{fi/T). These dissipative coefficients are given as a 
function of the temperature for five values of the baryon density: 0.2, 0.7, I, 
2 and 4 times the saturation density po. 

For the lattice eos, we made a fit of the shear viscosity from the results of 




(16) 



7 



LATTICE EOS 



0.3 



0.1 - 



0.2 - 












0.1 



0.2 0.3 



T [GeV] 



[31, 34, 35] which reproduces the characteristic behavior above Tc = 0.2 
GeV [36] and have a sharp increase around Tc [34]. (As a matter of fact, 
the cross section is expected to present a sharp change at a phase transition, 
due to the critical opalescence phenomenon. Then, the transport coefficients, 
which are in first approximation inversely proportional to the cross section, 
follow this behavior). 

2.4 Shock waves 

In the following, we will consider only the one-dimensional problem. In this 
case, the Navier-Stokes equations take the simple form 



M = {e + P)-f^v + K-f^{v^ + l)Vg-{4:ri + C)-f^v9/3 (20) 

with E = r°°, M = ^°^ = (7,0,0,7t;), 9 = dt{^) + d,{^v), Vg = 
vdt{fi/T) + d,{fi/T). 



dt{M) 
E 





(19) 
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Figure 6: Viscosity, 
thermal conductivity 
and K coefficient from 
Eqn. (4) function 
of temperature and 
baryon number [33] 
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We can see that in the simple 1-D case, the viscosity terms formally play 
the role of an effective pressure 

Pe// = P-{iv + 0^/3. (21) 

As 6 remains negative throughout the shock phase, the effective pressure 
gets higher than the thermodynamical one. We can also define an effective 
"velocity of sound" c^jj = Pejf/e. For too high values of this parameter 
- which can result from high values of the viscosities or too high values of 
the gradient -, the solution of the hydrodynamical equations would become 
problematic and ultimately would go out of the range of validity of the theory. 
It has been checked in our calculations that c^jj never exceeds I (It stays 
at the level of a 10% correction). The formal similarity between ideal fiuid 
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equations and one-dimensional viscous fluid equations is very probably at the 
origin of earlier interpretations of the experimental data, which needed in an 
ideal fluid model harder equations of state (with higher c^) than expected. 

In the conflguration of a heavy ion collision, a discontinuity in the ve- 
locities at the contact surface between two colliding nuclei results in the 
formation of shock waves which then propagate outwards. ([19] - [21]; [4, 5] 

)• 



Figure 7: Approxima- 
tion of the geometry 
of the nuclear colli- 
sion through two infl- 
nite tubes. The shock 
front (black bars) is 
moving with the ve- 
locity Vsh = C^maxl^eq- 



In Fig. 7 we visualize the basic features of the model. The participant 
region of the two colliding nuclei is approximated by two semi-inflnite tubes of 
matter. In front of the excited nuclear matter two shock waves are generated, 
moving apart and propagating into the incoming streams of projectile and 
target nuclear matter with the velocity ibc^/ueg- Here, Vpq is the velocity of 
the nuclei in the equal velocity frame. 

In the ideal case of two inflnitely extended nuclei which we flrst consider 
here, a stationary shock establishes. This should remain a good approxima- 
tion for large nuclei. In this case, the time derivatives vanish in the shock 
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frame, so that we obtain the set of equations 



p^v = 






--f'^vO = 


constant 




^^^0 = 


constant 


2 2 
= eo7r 



which after some rearrangement can be written as a set of differential equa- 
tions for the velocity and the ratio of chemical potential to temperature 

p^V = pQ^rVr (25) 

d,{plT) = [e^^lv.,{l-vv.,)-ev]l{^K) (26) 

£07>r[7'(l + ^^')(1 - VVr) - 1] + CT7'(c' " ^^') r 

= (4, + C)T'./3 

Here Sq and /Jq are the energy and baryon number density in cold nuclear 
matter, Vj. is the relative velocity of the inflowing matter and is related to 
the velocity in the equal velocity frame through 

Vr = [Veq + Vsh)/{1 + VeqVsh) (28) 

Vsh being the velocity of the shock. 



2.5 Asymptotic values 

For large enough nuclei, we can consider that the derivatives dz{v) and 
dz{p/T) vanish at z ^ oo. The density of the shocked matter reaches an 
asymptotic value e^^oo = £max- Its value its obtained by solving the system 
of equations (22 - 24) for v = Vsh^ £ = £max7 P = Pmax and the derivatives set 
to zero. In this way we obtain the velocity of the shock 

^sk = {ciJv,,) (29) 

with = P/e {c^ax — P™"^/e™"^). ^ This expression was derived for 
the case where the pressure in the unshocked matter is vanishing. If this is 

has to be distinguished from Cg = \/ dP/de\s=const (s:entropy-density) which is the 
speed of sound in the shocked matter. Although our P is not the speed of sound, we will 
use this loose way of speaking throughout the paper. 
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for some reason not the case, the right hand side of (23) and (24) have to 
be replaced by (eq + -Po)7r^r and (eq + Po)^lv1 + Pq respectively. Defining 



c. 



mm 



Po/co^ the shock velocity is then given as the solution of a 4 degree 
equation obtained by dividing (23) by (24), and replacing by its expression 
(28). From its solutions, we have to discard three as unphysical. We are left 
with the fourth solution 

V I, = — ■ — 1? -I- ■ 1? 

^ sn I max mm eq ' mm max eq 

+ ^U{l + ■ V"^ + (C^ ■ - C2 + ^2 _ . v"^ V] 

I w -Ly-i- I ^ratnj max eq ' \ mm max ' eq mm max eq J J 

/[2(l + cLJ^e,] (30) 



which is somewhat smaller than in the case ci, • = 0, and tends toward; 



s c 



irrespective the value of c^^^ when — > 1. The expression (30) we derived 
for the velocity of the shock is more general than the one derived in [18] 
{vsh = c^max) extcuds the model from ultrarelativistic applications to 
relativistic and semirelativistic problems ( SIS and AGS energies). 

The expression (29) for v^h also provides a lower limit in Eiah for the 
formation of a shock: 



V 



eq > ^raax =^ > m(l + c') / {I - c')' (31) 



e.g. for = 1/3, Eiah > rrin + 0.235GeV. Below this energy, we have a 
simple overlapp of the two densities rather than a shock wave. 

Behind the shock, the energy density reaches the asymptotic value 

£max = £obrM^ + ci^n)]/bl^sh{l + ciax)] (32) 

The value Cmax for the energy density behind the shock front is the same, 
whether there is dissipation or not. This behaviour comes from the assump- 
tion of infinitely extended nuclei. In that case, even for a large mean free 
path [i.e. a large width of the shock front), the instreaming matter would 
be stopped sooner or later. For more realistic nuclei with finite radii, this is 
evidently not the case anymore, (c./. §4) 

It is also interesting to express Cmax as a function of and the en- 
ergy available in the equal velocity system -E/afc, with (for c^^^ = 0) Veq = 



e.llq{vl + ^)l^ (33) 
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Maximal energy density reached behind the shocl< 



cl=0.2 
c^=0.6 



lattice (2) 
non linear QHD-1 




Figure 8: Energy den- 
sity reached behind 
the shock function 
of the laboratory en- 
ergy. The equation of 
state is here taken to 



One sees that the maximum reachable energy increases wit h Eiab and will 
be higher for softer equations of state, (see Fig. 8) 

Similarly, one has the asymptotic baryon number density 

Pmax = = P0^eq{Ve„ + C )/c (34) 

^shVsh 

From this, we also can derive the entropy produced across the shock 

{^maxi^ ~\~ ^max) Pmax Pmax) / Tjyiax (^^) 



3 Applications 

3.1 Width of the shock front 

The effect of the transport terms is to produce a sharp but continuous change 
of the variables across the shock front instead of a discontinuity. 

Let us first neglect the contribution of thermal conductivity. The accel- 
eration of the matter through the shock front is now given by 



dv 
dx 



{Ar, + 0/3 7'^; 



[{l + c^Yi^il-vVr)-!] 



(36) 
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and the energy density given from 

e = eoj'^Vril — vvr)/v (37) 

The characteristic length scale over which the change of the thermody- 
namical variables takes place - the width of the shock - can be evaluated 
through 

S = r {dvldx)-^dv (38) 

It is proportional to the viscosity, as can be read from Eqs. (36-38). 

The width 8 is plotted in Fig. 9. It is seen that it strongly depends on 
the eos, i.e. it increases with c^. In general we observe a decrease with the 
bombarding energy. 

The curve for the lattice model with phase transition presents an inter- 
esting peculiarity: After falling down to 1 fm, the curve rises again around 
-E'/a6=6.5 AGeV, corresponding to the appearance of the phase transition to- 
wards quark gluon matter. The rise is more spectacular when the viscosity 
also presents a sharp variation at Tc, but is still present if the viscosity is 
supposed to be constant. 

This behavior has two distinct origins: First, in the transition region, the 
density of the matter increases drastically with the collision energy. This 
corresponds to a quick variation in the parameter c, which shape is reflected 
in the curve for 8. Secondly the width increase at transition is enhanced by 
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dissipation. As a matter of fact, tlie viscosity is expected to increase strongly 
in this region, thereby increasing the width of the shock wave. 

Later we shall see that this behavior will influence the nuclear trans- 
parency in the phase transition region in a characteristic way. (see §4) 

3.2 Shock profile and the influence of the thermal con- 
ductivity 



Walecka 
c2=l/3 



> 

o 




Figure 10: Energy 
density and velocity 
profiles in the equal 
velocity system 



o'v® the density anfl veiociiy profiiKJ as seen in the equal velocity 



system for a bombarding «gy i?,.,=11.5 AGeV (typical at AGS) for three 
equally distant time steps. The plot was obtained by solving ( 36 - 37), i.e. 
by neglecting the contribution of the thermal conductivity. We can see that 
the maximal reachable energy density in the Walecka model is lower than in 
the case of a constant speed of sound 1/3. Also, the shock wave propagates 
faster for the QHD-1 eos. As a matter of fact, the average speed of sound 
in the QHD-1 model is higher than 1/3. Therefore (see 29, 32) e^ax is lower 
and is higher. 

Up to now, we have been neglecting the effects of the thermal conduc- 
tivity, as is widely done in hydrodynamics calculations with dissipation. We 
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would like to stress here that this assumption might not be justified. The 
numerical values of the transport coefficients are hardly known experimen- 
tally, however the microscopical models [3f, 33] rather predict the (shear) 
viscosity rj and the parameter K (related to the thermal conductivity, (c/. 6) 
to be of the same order of magnitude (see Fig. 6 ) 

Let us study the modification brought about by the taking into consid- 
eration of thermal conductivity. In order to do so, we rewrite (26 - 27) in a 
form similar to (36 - 37) in such a way as to separate the contribution of vis- 
cous dissipation, and the additional contribution of the thermal conductivity 
terms: 



e = eQ^lvr{l - v.Vr)lv - ^KVq/v (39) 

Since e increases across the shock and /u/T decreases with increasing e, 
the product KDq is negative. The additional contribution to e is then posi- 
tive. The additional contribution to dz{v) can be either positive or negative, 
depending on the sign of (u^ ~ c^)- Foi" high velocities it is positive so that 
we expect the velocity first to drop faster than without thermal conductivity. 
Then, when v has dropped below c, we expect the slope to decrease and v 
to return slowlier to Vgh. The global effect is an additional broadening of the 
shock front and a slight change in the shape of the profiles. 

In Fig. 11, we compare the shock profiles obtained with and without 
taking the thermal conductivity into consideration. The curves were obtained 
for a laboratory energy Eiah = 4 GeV, on the right for the non linear QHD- 
1 eos and on the left for the resonance gas eos. For both calculations, we 
used the parametrization of the viscosity and thermal conductivity given by 
Danielewicz [31]. We see that the energy density and velocity profiles have 
their width increased by say 15 to 20 %. The effect is much more remarkable 
on the temperature profile. The temperature has already increased much 
before a notable change has ocurred in the energy density, due to a better 
transport of T when the thermal conductivity is taken into consideration, 
whereas changes in T could only be brought about indirectly by the equation 
of state if only the viscosity was taken. 
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4 Dissipation and transparency 



As we already discussed the maximal density (Eq. 32) will only be reached 
for the collision of infinite tubes of nuclear matter, because the length of 
the colliding nuclei has to be larger than the width of the shock wave 6. 
Therefore the maximum density actually reached will depend on the size of 
the nuclei. 

If the width of the shock is larger than the diameter of the colliding 
nuclei the build up of pressure in the region behind the shock will stop in the 
moment when the shock is reaching the edge of the projectile or target. This 
gives a simple explanation for transparency from the viewpoint of dissipative 
fiuid dynamics. 



Shock profile at AGS energies (EOS Waleclca) 




Figure 11: 
baryon density profile 
in the case of a nu- 
cleus collision at 11.5 
AGeV. The horizontal 
lines are indicating the 
Lorentz contracted di- 
ameters of a Si and a 
Au projectile 



It also gives an estimate of the duration of the shock 

T ■ , (41) 

Veq + Vsh 

with R being the radius of the colliding nuclei. 

In Fig. 12 we compare the density profile of the shock front with the 
Lorentz contracted diameter of a Si and an Au nucleus. It is seen that in 
the case of the Si the maximum density is not reachable whereas in the Au 
case it is possible. One should also note the shape of the shock wave. Over 
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a long distance (~ 0.5/m) the density increases only very slowly and after 
that very sharply. This means that the projectile penetrates into the target a 
remarkable distance before something really happens. As already noticed, in 
the presence of a first order phase transition, the transparency effect should 
also lead to observable experimental consequences. The qualitative behavior 
of the transparency should follow the one of the width 6 (c.f. Fig. 9), 
i.e. we expect the appearence of a sharp increase around the transition 
temperature. However, due to the uncertainties in the determination of the 
transition density or temperature, it is not clear where this behavior really 
starts. Therefore one has to treat our prediction from ( Fig. 9) with care. The 
position of the bump will shift with the critical density. A closer investigation 
of transparencies in the range Eiab = 5-15 GeV could then give interesting 
clues as to the equation of state and for what critical temperature and density 
does the phase transition take place. ^ 



Duration of the shock 




Figure 12: Duration of 
the collision of two Au 
nuclei function of 
the laboratory energy 



The next two figures (14,15) display the energy density actually reached 
according to our dissipative shock model in a collision between two nuclei of 
mass number A in a fixed target experiment at energy Eiab- This parameter is 

^Let us stress again that the phase transition will however not actually occur, but 
rather be postponed to much higher laboratory energies due to the simultaneous increase 
of transparency 
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Figure 13: Energy 
den- 
sity actually reached 
in a collision func- 
tion of the mass num- 
ber of the nuclei A 
and the laboratory en- 
ergy Eiab- The mat- 
ter is described by 
the Hagedorn reso- 
nance gas model (see 
§3.1), the viscosity by 
Danielewicz's formula, 
A = 



defined from Fig. 16 as the value at x = and t = tmax = {RHr) I {veq + ^^s/i)- 
We displayed the two cases of a lattice eos and the resonance gas model. In 
the case of the non-linear QHD-1 model, the energy reaches almost its maxi- 
mal value eqn. 32 except for a narrow band A < 10, due to the irrealistically 
high value of the velocity of sound of this eos in the high temperature region. 
For the resonance gas model as well as for the lattice equation of state, we 
obtain similar shapes. There is a small structure around 3-6 GeV labora- 
tory energy, corresponding to the sharp rise and subsequent dip of the speed 
of sound in this region (see Figs. 2 and 5). The energy then rises slowly with 
A and Eiah for Eiah ~> 15 GeV or the lattice eos, the increase is sharper and 
starting from Eiah ~> 10 GeV for the resonance gas eos. As an example, 
for a Au-|-Au (resp. Si-|-Al) collision at AGS energy Eiah = 11-5 GeV, the 
energy density reached for our lattice eos is ~ 1.7 Gev. (resp. ~ 0.85 GeV). 
Finally we would like to point out that the values given here are probably 
overestimated, since (i) the definition taken here e{x = 0,t = tmax) is the 
peak value in Fig. 16 which is likely to be rounded off to smaller values in an 
actual hydrodynamical simulation and (ii) the plots displayed in Fig. 14, 15 
were made neglecting the thermal conductivity, which will give an additional 
contribution to the width of the shock front and thus further decrease the 
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actually reachable energy (see Fig. 11). 



Energy density actually reached (lattice eos) 



Figure 14: 
Fnergy density actu- 
ally reached in a col- 
lision function 
of the mass number 
of the nuclei A and 
the laboratory energy 
Eiab. The matter is 
described by the lat- 
tice equation of state 
(see §3.1), T] oc T"^, 
A = 
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5 Link with the full hydro dynamical simula- 
tions 

With the help of the parameter r defined in (41), we are able to show in Fig. 
16 what is to expect for the evolution of the energy density profile when the 
finite size of the colliding nuclei is taken into account. It is interesting to 
compare this to the result of a full hydrodynamical calculation (with rela- 
tivistic ideal fiuid hydrodynamics and a numerical viscosity of the order of 



The two models are in quite good agreement which each other; the only 
noticeable difference is the slope of the energy rise which comes from using 
different viscosities in the two models: the numerical viscosity used in the 
hydrodynamical calculation is an order of magnitude lower than that ex- 
pected in physical situations and used in our shock calculations. The shock 
model has the advantage of its simplicity, whereas the full hydrodynamical 
calculations are plagued with numerical problems originating in the size of 
the gradients involved at this stage. This argues in favour of using the shock 
model in order to modelize the early stage of a heavy ion collision. The results 



0.001 GeV.fm-2) [26]. 
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Comparison Dissipative Shock and HYLANDER (1-d) 




Figure 15: Density 
profile in the equal ve- 
locity system for finite 
nuclei (here, Au+Au) 
as given by our shock 
model (full line) and 
by a hydrodynamical 
model (dotted line) 
[30] 



of our shock calculations could then be used as an input for a subsequent 1-D 
(or possibly 3-D) hydrodynamical calculation in the expansion phase where 
the steepest gradients are already smoothed out and the equations better 
behaved. 

As an example, we have applied our results as input for a model with 
Bjorken-like expansion in the 1-D case or in the 3-D case with cylindrical 
symmetry. The 1-D hydrodynamical equation in the hypothesis of Bjorken 
scaling = x^/r, e = e(T), P = P{t), a(:= /u/T) = ^(t) is very simple: 



57+ — -V;5 = ° '«) 

Two remarkable properties of this equation should be noticed: (?) The ther- 
mal conductivity terms disappear from the equation for the 1-D Bjorken 
model (For a 2-D model with cylindrical symmetry and longitudinal Bjorken 
scaling, this is not the case any more) (m) Terms containing derivatives of 
the viscosities with r also disappear. 

If the equation of state can be put into the form P = c^e and c^, ?y, ( are 
assumed to be independant of r, (??) admits an analytical solution 

In this equation, Tq represents the instant at which the expansion phase 
starts and has to be equated with the duration of the shock phase (41). 
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e(To) is the maximal energy reached in the coUision (see Figs. 14-15). In 
Fig. 17, we display numerical results for the evolution of the energy density 
for Au+Au at Eiab = 4 GeV through the compression (shock) phase and 
subsequent Bjorken-expansion as given by (43). Four cases were considered: 
a small viscosity rj = 0.01 GeV/fm^c (full line), a viscosity in the order of 
magnitude of what is predicted by microscopical models i] = 0.05 GeV/fm^c 
(dashed line), a high viscosity rj = 0.1 GeV/fm^c (dotted line). These three 
cases show an important difference in the shock phase. Indeed, the finite 
width of the shock front arising from viscosity strongly modifies the rate at 
which the energy density is rising. On the contrary, in the Bjorken-expansion 
phase the rate of decrease of the energy density depends very weakly of the 
value of the viscosity, so that the continuous and dashed lines are practically 
undistinguishable on the figure. In order to make a comparison, we also show 
(dot-dashed line) the Bjorken expansion with an unrealistically high value of 
the viscosity rj = 0.5 GeV.fm^c, starting from the maximal energy density 
(i.e., for infinitely extended nuclei in this last case). The equation of state is 
assumed to be given by P = e/3. From this very simplified model, one would 
thus conclude that the infiuence of the dissipation is not very important in 
the expansion phase. However, the Bjorken scaling assumption is probably 
a too rough hypothesis. 
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6 Conclusion and outlook 



Let us first shortly summarize tlie results presented in this paper. 

We solved the 1-D relativistic dissipative hydrodynamical equations for 
a stationary shock. When only viscous dissipation is present, a shock width 
could be defined by an analytical expression. The width is decreasing rapidly 
with increasing laboratory energy; it is proportional to the viscosity and 
strongly dependent on the stiffness of the equation of state. It was shown 
to present a strong enhancement in the range 5 -12 GeV, where a phase 
transition to QGP would occur if the colliding nuclei were of infinite exten- 
sion, so that the matter behind the shock front had enough time to reach 
Cmax- However, the finite size of the nuclei prevents this maximal energy to 
be reached, instead we expect them to pass partly through each other. This 
could be an explanation for an unexpectedly high transparency observed at 
AGS energy (11.5 GeV) a few years ago [37]."^ 

The calculations were also made in case of a non vanishing thermal con- 
ductivity. It was found that the width of the energy density and velocity 
profiles are slightly increased. The temperature profile is much broader. The 
general conclusions as to the enhancement of the shock width and trans- 
parency are qualitatively not modified. 

Some questions which were not mentioned in the main text are worth 
being treated here. First, there is the problem of the stability of the shocks. 
Several non dissipative calculations ([22] - [25]) predict the shock to develop 
an instability and split instead into two shock fronts. In order to see if this is 
the case in our lattice eos, we display in Fig. 18 the Taub adiabat. The Taub 
adiabat is given by plotting the pressure P behind the shock as a function of 
a parameter x. When the system has a non vanishing baryon number /?, x is 
commonly defined as 



The equation of the Taub adiabat is obtained by eliminating the velocities 
from eqns. (22 - 24) with the gradients set to zero. The equation of the 

^According to Fig. 12, we would estimate the transparency to be K ~ p(Si) / p(Au) = 
0.66. This is comparable with the result of [37] who find that 11 out of 28 nucleons passed 
through the target without being affected by the collision, corresponding to K ~ 0.61. We 
do not think however that this coincidence should be taken too seriously. 




(44) 
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Taub adiabat is then found to be (the index 1 (2) describes the unshocked 
(shocked) matter respectively) 

(£2 + P2)x2 - (ei + Pi)xi = (P2 - Pi){x2 + Xi) (45) 

The quantities before the shock (index 1) and the equation of state behind 
the shock being known, (45) is an imphcit equation for P2{x2). 

However, when the baryon number is vanishing such as is the case for 
our lattice eos, x has to be defined in an other way. It was proposed by 
Danielewicz and Ruuskanen [32] to use instead 

x = [{e2 + P2hM]/[{e^ + P^hy,] (46) 

One obtains in this case formally the same equation as in (45) with the 
replacement Xi — )• 1. 

When moreover the pressure in the unshocked region vanishes, as is the 
case for cold matter at nuclear density (and as was assumed throughout this 
paper), the Taub adiabat formula simplifies to a transparent form 

P = 

-'■ max / I \ 

AXjXQ) 

max 



(47) 
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with Xq = £o/pI for nonvanishing baryon number, 1 otherwise. 

In order for the system to be stable against the formation of a double 
shock, the pressure has to be a monotonously decreasing function of x be- 
tween X^ax = Xo (Vs/Vr) = Xq {c^ / V^^) {v^^ + C^) / {1 + C^) and Xq. OuC SCCS 

that this is the case for the non-linear QHD, the Hagedorn resonance gas 
model and the lattice model for parameter set 2. For the lattice eos with 
parameter set 1, the pressure is not a monotonous function of x anymore , so 
that a double shock could form in principle. However, this is not so simple in 
the case of a dissipative shock because it smears out the sharp separation in 
two shock fronts seen in the case of an ideal fluid [22]-[25] , so that what we 
observe is a single shock front with an enhanced width rather than a double 
shock. 

We would also like to say a word about the principle question raised by 
e.g. Olson and Hiscock [16], (see also [12] - [15]): It has been argued that 
the relativistic Navier-Stokes equations in their usual form (eqs. 5-8) suffer 
from acausality and instability, ("acausal propagation of heat"). We are 
conscious of this difficulty. However, we think that, for practical purposes, 
the relativistic Navier-Stokes equation still has a long life expectancy. Of 
the alternative theories proposed, the method of moments with inclusion of 
higher order derivative terms in the right hand side of (eqs. 3,5,6) d la Israel 
and Stewart [13] or non local terms [14], none of them is yet 100% satisfactory, 
and would need further investigation. They have the big disadvantage of 
introducing additional complexity through the higher order derivatives, while 
even numerical hydrodynamics using the simpler Navier-Stokes equations 
are still in their infancy. 9 additional parameters have to be introduced 
in the Israel theory as coefficients of the higher order terms, which values 
are not or very bad known theoretically nor measurable experimentally in 
a foreseeable future. The uncertainties on these parameters again would 
introduce uncertainty in the theory. 

There are other consequences of the finite size of the nuclei which we did 
not mention in the main part of this paper. The two most important ones 
are (i) The evaporation of particles in a preequilibrium stage through the 
surface [38] (ii) The reduction of the transport coefficients from the finite 
size effects (see e.g. [39]). As a matter of fact, the transport coefficients 
displayed in Fig. 6 were calculated in a model of infinite nuclear matter. In 
finite matter, the particles situated near the surface of the nuclei have their 
mean free path limited in one direction and collide with the potential wall of 
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the nucleus rather than with other particles. This results is smaller values 
of the dissipative coefficients (as measured for example from fission necks or 
"hot spots") near the surface, i.e. (see e.g. Fig. 11) in a region of smaller 
temperature, thus counteracting the rise at low T (?y oc l/T^) exhibited by 
the kinetic solution of [31, 33] (see Fig. 6). 

Finally, we would like to mention other important consequence of the 
model concerning the emission of Bremsstrahlung from the acceleration through 
the shock. From the characteristics of the Bremsstrahlung spectrum, one can 
have access to the energy density reached, the width of the shock front, and 
the duration of the shock phase. This in turn may give information about 
the stiffness of the equation of state and dynamical properties of the matter. 
This question is currently under investigation and will be the subject of a 
forthcoming paper. 
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Figure 18: Energy density, velocity and temperature profiles in the rest frame 
of the shock, for the non linear QHD-1 and the resonance gas equations of 
state, at Eiab = 4 GeV. The profiles obtained by taking the both the viscosity 
and thermal conductivity into account (full line) are compared with those 
obtained by considering the viscosity only (dashed line) 
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